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Abstract 

We study geometry of the phase space for finite dimensional dynamical systems 
with degenerate Lagrangians. The Lagrangian and Hamiltonian constraint formalisms 
are treated as the different local-coordinate pictures of the same invariant procedure. 
The invariant description is given in terms of geometrical objects associated with the 
structure of foliation on the phase space. 
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I. INTRODUCTION 



Since the appearance of the fundamental work of Dirac fll], a certain progress has been 



made in understanding geometry which underlies constrained dynamics || fll3[ . For a 
system with regular (nondegenerate) Lagrangian an appropriate phase space in the case of 
Lagrangian (Hamiltonian) formalism is the tangent bundle TQ (cotangent bundle T*Q) of 
the configurational space Q. In consideration of degenerate Lagrangians it is traditionally 
assumed that the phase space of the system is a subspace of TQ (Lagrangian formalism) 
or T*Q (Hamiltonian formalism). Tangent and cotangent bundles are treated separatively 
within the framework of the constrained dynamics, being related by the standard Legendre 
transformation. In this paper another point of view is advocated. 

Given a Lagrangian I 6 C°°(TQ), the action S can be written as a functional 



S[a] = I *V 



acting on the sections s of the bundle TQ(g)R — > R, where X e R is a compact set, 
ip = —Edt + 81, Oi = ^dq is Cartan - Liouville one-form and E is the energy of the system. 
The Legendre transformation L: TQ — > T*Q is a diffeomorphism for nondegenerate I, so 
the transition to the Hamiltonian formalism can be viewed as a change of variables in (1), s 
being the section of T*Q (g)_R — >• R, <p — —hdt + 9, where 9 = pdq, Oi = L*9 and h is the 
Hamiltonian, E = L*h. 

For degenerate Lagrangians the map L is in general neither surjection nor injection. In 
this case the pleasant symmetry between the Lagrangian and Hamiltonian formalisms is 
broken, since the image of TQ by L now is only a subspace of the cotangent bundle. If 
one wishes to proceed in the usual way, the machinery of Dirac - Bergmann algorithm [p], 
is applicable either in local-coordinate approach or by using the invariant (coordinate-free) 
description in the spirit of Gotay - Nester's work [§, |, 

The invariant description of the constrained Lagrangian dynamics begins with the basic 



notion of foliation on the tangent bundle and should in principle end up with the well-defined 
procedure of the reduction of the system to the submanifold of the lower dimension (which 
itself must admit a tangent bundle structure) with no second class constraints left in the 
theory. The next step is to compare the result with that one of the reduction procedure in 
the Hamiltonian formalism and check that the reduced submanifolds are connected by the 
nondegenerate Legendre transformation. For the specific choice of the Lagrangian function 
the reduction procedure of that kind has been described in |2]J by making use of the tangent 
bundle geometry formalism developed earlier in pO . 



The Hamiltonian reduction is controlled by the standard symplectic structure on T*Q 
given by the closed nondegenerate two- form uj = dp A dq, whereas the analog of u in the 
Lagrangian formalism, ui = L*uj, is degenerate for degenerate £. 

The usual treatment of Lagrangian dynamics starts from the invariant form of the Euler 
- Lagrange equations on TQ, 

f = X E f, (2) 
/ G C°°(TQ), where for any g G C°°(TQ) X g is defined by 

ix g ui = -dg. (3) 

One may note that the relation (3) also plays the key role in geometric quantization [ITS], ITS]. 



In this paper we pay more attention to the action (1) itself. In the case of degenerate 
Lagrangians the bundle TQ ® R —>■ R appears not to be an appropriate phase space and 
should be replaced by the bundle W ® R — > R, where W = TQ T*Q is a Whitney sum of 
TQ and T*QQ . The equations which come from the variational principle can still be written 
in the form of (2) but the relation (3) must be deformed taking into account the structure of 
the foliation on the original phase space W. This modification allows us to include into the 

scheme those functions (such as the second class constraints) which are not constant along 

^^After this work was completed, we learned that similar approach was previously considered (in local 
coordinates) by W. Kundt [E5J and by R. Skinner and R. Rusk (using coordinate-free formulation) p6| 



the leaves of foliation. Then the equivalence of the Lagrangian and Hamiltonian formalisms 
becomes apparent and the constraints in both pictures are in one-to-one correspondence with 
each other. For nondegenerate t the original phase space W decouples trivially into TQ and 
T*Q at the first step of the constraint algorithm procedure. 

It follows that the most general situation one encounters in the invariant description of 
the classical dynamics is that of the presymplectic manifold rather than the symplectic one. 

The organization of the paper is as follows. In Sec. II a brief review of the invariant 
formulation of the Lagrangian and Hamiltonian dynamics is given for the systems with 
nondegenerate Lagrangians. The degenerate case is treated in local coordinates in Sec. Ill 
and the invariant formulation is provided in Sec. IV. Several illustrative examples are given 
in Sec.V. The Appendix contains definitions of some geometric objects used. 

II. NONDEGENERATE LAGRANGIANS 

A. LAGRANGIAN FORMALISM 

For a given Lagrangian l(q, q) the action 

S[q(t)}= Je(q,q)dt (4) 

with the appropriate boundary conditions can be deliberately changed to 

S[qM = fdt(e(q,v)-v i ^ + ?^Jdt (5) 

to be considered as a functional on the tangent bundle TQ of the configurational space 
(smooth manifold) Q. We denote the independent local coordinates on TQ as (q l ,v l ), i = 
1, ...N, N =dim Q, l(q, v) G C°°(TQ). 

The Lagrangian I is called nondegenerate (regular) if the matrix T^- = d 2 l/dvidvj has 
maximal rank R = N. Otherwise the Lagrangian is called degenerate (singular). Note that 
R does not depend on the particular choice of the local coordinates. 
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The variational principle applied to (5) gives 

6S_ _ dl_ __ d 2 l k d 2 £ . k _d_ ( dl_ \ 

6qi - dqi - Q v k dq i v + Ov k dq iq ~ Jt \drf ) ~ ' ^ 

§ = T r 0-v>) = O. (7) 

For nondegenerate / the last equation implies q 1 = v l . Then (7) gives the Euler - Lagrange 
equations. One can define a special vector field X E e X (TQ) corresponding to (6), 



where 

M_ __ dH k 
013 dvidq kV {) 

and T tk T k j = 5^ such that for any / e C°°(TQ) its time derivative is given by 

/ = X E f. (10) 

In particular, for energy E = — i we have 

E = X E E = 0. (11) 

The invariant description can be given as follows. Let Q be a finite-dimensional smooth 
manifold, (TQ,tt,Q) - its tangent bundle, i G C°°(TQ). The Cartan - Liouville one-form 
6 L G A X (TM) is defined by 

6 L = d£o r, (12) 

where r is a canonical type (1,1) tensor field on TQ j|, ^3|. Then the closed nondegenerate 
(for nondegenerate t) two-form u L = d6 L provides TQ with the symplectic structure. In 
local coordinates 

e L = ^~dq\ (13) 



u L = Tijdtf A dq j + -Mijdct A dq>, (14) 




For /, g e C°°(TQ) one can define the Lagrangian bracket by 



{f,g} L = u L (X f ,X g ), 



(15) 



where the vector field Xf E X (TM) is defined by 



ijr^L = -df. 



(16) 



Since cj l is nondegenerate, the map f2 x (T(5) — > A" (TQ) (16) is an isomorphism. In local 
coordinates we have 



where T' lk Ykj = Sp M lk = P fe M^r' J . This bracket is antysimmetric. It obeys Jacobi identity 
since du>L = 0. The energy E can be introduced as a smooth function on TQ by 



where Z = {X £ X(T Q):tX = A}, A is the Liouville vector field. 

Putting all this together, we can write the action (5) in the invariant form (1). For 
the energy (18) one can define the corresponding vector field using (16). A simple local- 
coordinate calculation shows that this vector field is equivalent to the field Xe (8) which 
originates from the variational principle. Therefore for any / G C°°(TQ) we can write 




(17) 



e = i z e L - £, 



(18) 



f = X E f = df(X E ) = u L (X E , X f ) = {E, f} L . 



(19) 



In particular, 



E = dE(X E ) = 0. 



(20) 
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B. Hamiltonian formalism 

The transition to the Hamiltonian form of dynamics is provided by the Legendre transfor- 
mation L : T q Q — > T*Q 

d£ 

(q,v) -> (q,p), pi = —, (21) 

(q,p) G T*Q which is a diffeomorphism for nondegenerate /. In the new variables the action 
(1) becomes 

S[q,v] -> S[q,p], 

S = I s>, (22) 

where ip = -hdt + 6, ip G tt\T*Q®R), X C i? 1 , %,p) = (L*) _1 (£?) is the Hamiltonian, 
^ = pj(i<f is the canonical one- form on T*Q, 9 L = L*(9). 

The variational principle applied to (22) gives the standard Hamilton's equations which 
can be associated with the vector field X h E X (T*Q)(the counterpart of X E (8)), 

X h = — — - —— (23) 
dp 1 dq l dq l dp 1 

The field X h can be defined invariantly by 

ix h oo = -dh, (24) 

where u = d9 is the nondegenerate two-form which provides T*Q with the standard sym- 
plectic structure. Then V/ G C°°(T*M) we have 

f = X H f = df(X H ) = cu (X H , X f ) = {H, /}, (25) 

where {/, g) is the Poisson bracket. In particular, 

h = dh(X h ) = 0. (26) 
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III. DEGENERATE LAGRANGIANS 



For degenerate Lagrangians the Legendre transformation (21) is not a bijection. This means 
in particular that R= rank T^- < iV so the condition x % = v % is valid now only for a part of 
coordinates, i.e. for x a = v a , a,—l,...R. The transition to the Hamiltonian formulation can 
no longer be considered as a change of variables in the action (1). To cure this situation one 
should extend the phase space by adding N — R new independent variables p M (/z=l,...iV — R) 
and write the action S as follows fL5 



r ( dl d£ \ 

S[q,v,p,} = J dtU(q,v)-v«— + q a —+ P ^-v»)j (27) 

which can also be presented in the form 

S=[(pi, (28) 



with ip i = —Eidt + di, where 

dl 

E l = v"—+p^-l, (29) 



dv° 
dl 

e^—dq^p^. (30) 



The variation of S gives the equations 



3/ 



dl d' 

v a = r afe ^ - r a6 ^-^ - r*r^, (3i) 

dl 

i» = W . (32) 

q i -v i = 0, (33) 

*-£ = °- < 34 > 

The last equation is a constraint in the space of variables q l , v i ,p fi . The time derivative of a 
smooth function f(q l , v\p^) can be written now as 

f = Y^f + v^f, (35) 



where 



dl 

^ = *V - = 0> ( 36 ) 



<9g* 6 <9t> a ~^ dq^ dp^ ' ^ ^ 
dt d 2 £ 

ab = W~ W'' (38) 

n w = ^ < 39 > 



and 6^ are the null- vectors of T^-. 



The action (28), however, does not have an invariant form. For the sake of invariance 
one should write S as 

S[q,v,p] = J dt(l-p t v l + p t q l ). (40) 
Application of the variational principle gives dynamical equations 

P. = |, (41) 
q< = v i (42) 



as well as the primary constraints 



dl 

&= ft -^=0. (43) 



The action (40) can be considered as the basic one for both the Lagrangian and the Hamil- 
tonian formulations for either degenerate or nondegenerate Lagrangians. 

Equations (41)-(42) imply that the evolution of any smooth function f(q,v,p) is given 

by 

f = Yf + v%f, (44) 
dl 

h=P*-W = °> (45) 



where 



d dl d 

Y = v l — + __ (46) 
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r > = w (47) 

Since Yi<pj=-Tij, the Dirac's compatibility condition fa = leads to determination of R 
accelerations 

v a = (Yfa)T ba - v tl T llb T ba (48) 

as well as to the new constraints 

*f = % (Yfa) = 0. (49) 

At the same time a part of constraints in (45), namely those for which 

det(Y a fa) ^ 0, (50) 

Y a = can be resolved. Because of the condition (50), we are free to choose either v a or p a 
as the local coordinates on the resulting surface M x = {q\p i ,v' t : fa = 0,a = 1, ...R}. As a 
result of this choice one gets either Lagrangian or Hamiltonian formulation, correspondingly. 

A. LAGRANGIAN FORMALISM 

In the Lagrangian scheme the system (44)- (45) on Mi becomes 

/ = yw/ + ww/, (5i) 

0Jm x = 0, (52) 
\ Ml = 0, (53) 



where and F^ 1 ^ have been defined in (37)-(39). The action S (40) confined to M 1 is 
given by (27). We will call the primary constraints for which (50) (and consequently (48)) 
holds the primary constraints of the second class. The rest of them are called the primary 
constraints of the first class. 

Next the condition fa^ = is to be verified. This gives R\ determined accelerations, 

v« = -(YW<$)V -&»{YV<$)y , (54) 
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and a new set of constraints, 

« = 6g" = 0. (55) 

Here 7^ = Yjpffl, R l= mnk ai =l,...R u m=l,...N - R - R l and bf/ = - 

^ai7 ai6l 7biMi are ^ ne nu ^ vectors of 7^. Again, among the secondary constraints (53) there 
are R\ constraints of the second class (det7 a 5 ^ 0) as well as iV — R — R x those of the first 
class. The system therefore can be reduced to the surface M 2 C M 1 of the second class 
secondary constraints 0^ = (the corresponding number of the primary constraints (Ri) 
must also be resolved, since ai ~ 4>*£})- The evolution then is given by 

f = YWf + &"YWf, (56) 



|m 2 = 0, (57) 
<f> { £ \m 2 = 0, (58) 
4? \m,= 0, (59) 



where 



\M 2 : 



(60) 



= im 2 • (6i) 

The following iterations are straightforward. After the k-th step the system becomes 

/ = y(*)/ + »«-'yW i / ) (62) 

plus the set of constraints 

^L, k= °» ( 63 ) 

where a = 0, 1, ...A;, 0^ = 

= (y(*-D _ (Y^M k k :?) (fc 7 1)6 " ia " 1 F^) k, (64) 
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M k = {xe M fe _! : (p^ix) =0,a = 0, ...k}, M k C M fc _x C M k ^ C---CW. 
The iteration process terminates, if: 

1) for certain k we have 

Rk = N — Ri — ■ ■ ■ — Rk-i, 

(k) 

where R k = rank 7 fe fe _ 1 a fe _ 1 - in this case all the accelerations are determined and the dynamics 
confined to the final constraint surface M k is totally fixed. 

2) all the constraints of generation k + 1 are reducible to those of the previous generations 
or are identically zero. In this case N — R\ — ■ ■ ■ — R k accelerations remain undetermined 
and the system posesses certain "gauge freedom" . 

B. HAMILTONIAN FORMALISM 

Let us rewrite the system (44) - (45) in a slightly different form 

f = {EJ} + v i Y i f, (65) 

<t>i = 0, (66) 

where E(p, q, v) = PiV 1 — I and {/, g} is the usual Poisson bracket generated by the two-form 
dp A dq (which is of course degenerate in the space of variables p,q,v). As it was in (54) - 
(55), we determine R accelerations, 

t; a = {E,0 fe }r 6a -^r^r ba , (67) 

and a set of new constraints, 

4>f = 6j, {E, = 0. (68) 

The condition (50) allows us to resolve a certain part of constraints in (66). To produce the 
Hamiltonian scheme, one should choose p a rather than v a as a set of local coordinates on 
the surface M 1 of the resolved constraints. Due to the identity = T^ a T ab T bu the Routh's 
function 

R(q\p a ,vn = (l-PaV a ) \ Ml (69) 
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is linear in 

R=-h(q\ Pa )+v»M<l\Pa) (70) 

(here ^>^(q\p a ) = dl/dv^ \m 1 ), so the action (40) being confined to Mi in coordinates 
p a , q\ reads 

S{q\Pa,v»] = I dt(-h T + Pl q l ), (71) 

where 

h T = h(q\p a )+vy fl , (72) 

Correspondingly, the system (65) - (66) becomes 

/ = {K /} + /} + t>%/, (73) 

M |m x = 0, (74) 
k= 0, (75) 



Note that (28) and (71) are equivalent. The transition between the Lagrangian and the 
Hamiltonian formulations is a smooth change of variables in the action, as it was for non- 
degenerate Lagrangians. There is no need to add voluntarily the primary constraints M to 
the function h(q\p a ), since the old Dirac's construction hr naturally appears in (71) and in 
(73). 

The condition 0^=0 leads to the determination of some accelerations, 

t/ 11 = - {h T , 0^} 7 6iai - v^ lbl l biai , (76) 
and to the new constraints, 

«=&?Ke ) }=°- ( 77 ) 

Here 7^=9(0^ \M 1 )/dv ,i - Now we can reduce the system to the surface M 2 C M 1 of the 
second class secondary constraints (/>$ = 0. In the Hamiltonian scheme we can do this 
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explicitly by expressing v ai : 

v a ^ = - {h, <f> bl } 7 biai - v" 1 , 6l } 7 6iai • 
Then the system (73)- (75) on M 2 can be written as 

|m 2 — 0, 

< } Im 2 = o, 

where 

{f,g}l = {f,g} - {f,<f> bl h biai {K,9} 

is the Dirac bracket in respect to the second class primary constraints. Next iteration 
us another set of the determined accelerations, 

v a > = - {Mg>}* tf™ ~v vi {<t> vl A { S}\ 7 )62<12 -v" 2 7^ 2 l b2a \ 

velocities, 
where 



and constraints 

The system (79)- (82) then becomes 

f = {h,f}l + v"> {<!>„, m + vVYnf, 

4>tX 2 \m 3 — 0, 
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1^=0,0 = 1,2,3, (90) 



where 



i * \ ) , , , 4 r ,11* W 



H/^K-f/,^}! V {0a 2 ^} 1 + {^,0£ ) } 1 V {0a 2 ,/K- (91) 

The iterative procedure continues until either all acceleration ( and velocities) are determined 
(and the final phase space is a submanifold of T*Q) or the constraints of the k-th generation 
are reducible to the previous ones (or are identically zero). 

IV. INVARIANT FORMULATION 

Let M be a real smooth N- dimensional manifold, (TM, 7T\, M) its tangent bundle, (T*M, 7r 2 , M) 
its cotangent bundle. Let W=TM $T*M denotes a Whitney sum of these two bundles 
(Fig.l). 

W 

Pi / \ P2 

TM T*M 

TTl \ / ^2 

M 

Fig.l 

Consider p\Q e Q l (W), where 9 is the canonical one-form on T*M. To simplify notations 
we denote p\Q by 6 again. The two-form uj = d6 is closed and degenerate on W, Ker 
u = X (TM)c X(W). The elements of Ker uj form an integrable distribution T> on W thus 
defining a foliation F of W, the leaves of F being the maximal integral manifolds of T>. The 
foliation F has codimension q = 3N— rankD = dim Coker u=2N. We assume that the 
space of leaves C has a structure of manifold, so F is a regular foliation [JD]] of codimension 
g. Now consider the de Rham complex of F. We have Si L W/£=T*W/£, WW / £=A l tt l W / £, 
nW/C=®i> &W/C. The first order differential operator d g : fiW/£ -> fiW/£ has the 
following properties |22| : 

i) d, : C^W) -> fi 1 !^//: is a composition C°°(W) S WW -> T*W/£, 
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ii) d q (u p A cu s )= d q cu p Acu s + (-l) p uj p A d q u s , uj p E tt p W/£, 

iii) d\ =0. 

Note also that do = d : VtW — * QW. In any Frobenius neighborhood U C W one can 
choose the local coordinate system (x 1 , ...x d , y d+l ...y d+q ) such that T>u is generated by the 
set {d/dx 11 }, /i=l,...d=rank V. Then the action of d q on any a E il p W/C, 

a = a fll ...^ lp (x,y)dx 1 A • • • A dx» p , 

can be described as follows 

= da K..4»( x >y) dx u Adx i A ... A dx » P 

ax v 

/i, ...d. In other words, provided {Ki} d =1 is a basis in V and {#*}^ =1 is a dual basis in 
£>*, the action of d q is given by 

d q a = [Kia h ... ip ) 6 i A 1 A • • • A 6 ip . 

Let us denote ujj = d q f, f E C°°{W). Now we are ready to proceed with the invariant 
formulation of dynamics. 

Let Z = {X E X(W) : r(puX) = A}, where r is the vertical endomorphism, A is the 
Liouville field. Then for any Lagrangian I E C°°(TM) the energy E E C°°(W) is defined as 
follows 

E = i z 9- p\l. (92) 
The action is a functional on the smooth sections of W R— > R , 

S[s] = [ 8 *<p, (93) 

where X E R is a compact set, (p = —Edt + 9, <p E ^(W <g> R). Since W is only a 
presymplectic manifold, the relation (16) between one-forms and vector fields on W now 
makes sense only for functions constant along the leaves of foliation F. In particular, (16) 
cannot be used for the second class constraints of any generation. However, in the local 
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coordinate approach to the constrained dynamics we do not encounter such restriction. To 
cure this, we deform the correspondence (16) and put it to be 

t Xf uj = -df + u q f , (94) 

where q is the codimension of the foliation generated by Ker uj. Then the system of equations 
(44) - (45) which follows from the variation of S can be written in the form 

f = df(X E ), (95) 

w q E = 0, (96) 
where X E is defined by (94). Using (94) we can rewrite (95) as 

f = u(X E ,X f )+u q f (X E ). (97) 

Note that (94) -(95) give E = 0. One can decompose X E as 

X E = Y + K, (98) 

where K E V is an arbitrary vector. Locally, K = a l Ki, where {Ki}f =1 is a basis in V. 
Since V E X V (TW), the yet undetermined multipliers a % are interpreted as accelerations v l 
in any given chart. The equation (96) then reads 

u% = (KiE) d { = 0, (99) 

with 8 l (Ki) = 5^. This defines a subspace S C W by 

(f>i = dE (Ki) = 0. (100) 

For compatibility of (95) - (96) one must have therefore X E E TS, i.e. 

dfa (X E ) = 0. (101) 

If S is transversal to the leaves of F (in this case V n TS = and rank | d<pi(Kj) \=N), the 
vector field X E is completely fixed by the condition (101) in the sense that all coefficients 
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a % are determined. Then we can confine the system to S and choose either p's or v's as the 
set of local coordinates on S thus generating either Hamiltonian (S = T*Q) or Lagrangian 
(S = TQ) form of dynamics. 

In the generic case, however, VnTS=Vi ^ 0, rank | d(j)i(Kj) \— Ri < N, so only R\ out 
of N accelerations are determined by (101). Instead of S we can choose now an intermediate 
transverse subspace M 1 (S C M 1 C W) such that d 1 =dim (V U TM X ) = N - R x . The 
two- form uji = uj \m 1 is degenerate, dim Ker lo x — d\- Now the scheme repeats itself. The 
system (95) - (96) becomes 

f = df(X El ), (102) 

c4 k= o, (103) 

< \m,= 0, (104) 

where E\ = E \m 1) qi = is a codimension of the foliation of Mi generated by 
This system can be expressed equivalently in terms of either Lagrangian ( (51) - (53) ) or 
Hamiltonian ( (73) - (75) ) local coordinates. A new subspace S\ C Mi defined by (103) 
arises here. The requirement X El C TS\ gives R2 < N — R\ determined accelerations. We 
choose M 2 such that d 2 = dim (V n TM 2 ) = N — Ri — R 2 . Now uo 2 — uj\ \m 2 is degenerate 
and dim Ker u 2 = d 2 . The dynamics is governed by the system 

f = df (X E2 ) = u 2 (X E2 ,X f ) + uf (X E2 ) , (105) 

u% \m 2 = 0, (106) 

< \m 2 = 0, (107) 

< \m 2 = 0, (108) 

The sequence stops when for certain n S n is not foliated, i.e. when P n _i n TM n _i = and 
<jj n is nondegenerate. In this case the dynamics on S n = M n C W is fixed and one has 

f = u n (X En ,X f ) (109) 
17 



for any / G C°°(S n ). As an alternative, one may discover, however, that for some n the 
condition Xe ji _ 1 G TS n -i is automatically satisfied while S n -i is still foliated. This means 
that N — Ri — . . . — Rn-i accelerations remain undetermined and the evolution is described 
by the equation 

f = u n (X En ,X f )+ujf(X En ), (110) 
together with the final set of the first-class constraints 

"!k>=o,...< U„=o. (in) 

The invariant geometrical treatment of the first-class constraints case including gauge-fixing 
procedure was given by L.D. Faddeev ||. 

V. EXAMPLES 

Example 1: Let us consider the simplest nondegenerate Lagrangian, 

l = Y'U(q). (112) 

We start with the phase space W = TR 1 ($T*R l . The canonical two-form u = dp A dq is 
degenerate on W, Ker to— span{<9/<9t> }. The energy E G C°°(W) (92) is given by 

E = pv-^ + U(q). (113) 

The primary constraint one- form uje can be written as 

dE 

uje = d„E = -r—dv = 6dv = 0, (114) 
ov 

where (j) = p — v. The relation (94) gives the evolution vector field X E on W, 

d dU 8 . d , iic , 
X E = v— - — — + v—. 115 
oq oq op ov 

Thus the dynamics on W is given by {/ = X E f, 4> = 0}. The compatibility condition (101) 
allows to determine the acceleration v = dU /dq. This means that <fi is the primary constraint 
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of the second class. Therefore, we reduce the system to Si = {x G W : 4> = 0}. One can 
choose either (q,v) or (q,p) as the local coordinates on Si. In the first case one gets the 
Lagrangian scheme, 



v 2 



Ei = E\ Sl =- + U(q), (116) 

X - V ^--^IL1- (117) 

El dq dq dv J 



in the second case - the Hamiltonian description, 



2 

h=E\ s = ^ + U(q), (118) 

9 dU d . . 

X h = p—-— — . 119 
oq oq op 

In this case W trivially decouples into TQ and T*Q. 
Example 2: Consider the Lagrangian (Ref. [|14]1 ) 

l = V \-v 2 qz- (120) 

The initial phase space is W=TE? 0T*i? 3 with local coordinates qi,Pi,Vi, % = 1,2,3. The 
fundamental two-form u = dpiAdqi is degenerate on W, Ker u = span {d/dvi, d/dv 2 , d/dv 3 }. 
The energy E G C°°(W), 

E = piVi +P2V2 +P3V3 ~ ]^vl + v 2 q 3 , (121) 

defines (by (94)) the evolution vector field X E = Y + K, where 

Y = v~ - «a A (122) 
oqi op 3 

K ="i (123) 

with yet undetermined multipliers Vi. The one-form (96) of the primary constraints is given 
by 



dE 

u) E = d q E = -T^dvi = (pidvi = 0, (124) 
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where <pi — Pi — v 1, 02 = P2 + <?3, 03 = Pa- The compatibility condition (101) determines one 
of the accelerations, i)\ — 0, and produces two new constraints, 

# =^s = 0, 3 1} =v 2 = 0. (125) 

Thus the primary constraint 0i is of the second class and can be resolved. On the interme- 
diate transverse subspace M x = {x G W : 0i = 0} the evolution is given by 

; df df df df , . 

J =v i* V ^ VV2 7> ht; 3^— , (126) 

oqi dp 3 ov 2 dv 3 

together with the set of constraints 

M 2 = {0 2 = 0, 03 = 0, 0^ = 0, 0« = 0} . (127) 

Since = ^3 = 0, = V2 — 0, all the secondary constraints are of the second class and 
the final constrained submanifold S2 = M2. For any / G C 00 ^) one has then / = Vx^-. 
The dynamics is totally fixed. 

To consider the Hamiltonian form of dynamics, one should choose p\ instead of v\ as the 
local coordinate on M\. Note that 

E \hh= + ^2<?3 + P2V2 + P3V3 = h + v 2 <p2 + v 3 (p 3 . (128) 
The equation (126) can be written in the form 

/ = {K /} + V2 {02, /} + ^3 {03, /} + V2^~ + (129) 

OV2 dv 3 

together with the set M 2 (127). The compatibility condition again gives i>2 = 0, v 3 = 0. 
One can make the canonocal transformation, Pi = pi, Qi = q±, P 2 = P2, Q2 = Q2 + P3, 
P3 — 03 = Psi Q3 — Q3+P2 = 02- In terms of these new variables S 2 = {P2 = 0, Q 2 = 0, P 3 = 
0, Q 3 = 0}. One has h | & = P 2 and / = {h, /} for any / G C°°(S 2 ). 

Example 3: Let / = \ (qiv 2 — q 2 V\ — q\ — qV) (Ref. |2^| ). The two-form u = dpi A dqi, 
i — 1, 2, is degenerate on W = TP 2 0T*P 2 , Ker oj = span {d/dvi, d/dv^}- The energy is 

P = Pif 1 + p 2 f 2 - I (130) 
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and the evolution vector field Xe is given by 

d d fv 2 \ d fvi \ d d d 

X E =V 1 ~ h V 2 - Ht-?! -"92 7 1" «1 a V ^277—, (131) 



dqi dq2 V 2 / dpi V 2 / dp2 dv\ dv 2 
The set of primary constraints is defined by 

uje = d q E = (fiidvi = 0, (132) 

where 0i = Pi + 92/2, 02 — P2 — 9i/2. The condition (101) gives only the new constraints, 
01 = v 2 — q± = 0, <f>2 =V\ + q2 = 0. The next iteration (J>i = O,0 2 ^ — 0), however, allows 
to determine i)\ = —v 2 and V2 = V\. Therefore, the final constrained submanifold is 

S = {0 1 = O,0 2 = 0,0^ = 0, 0^ = 0}. 

The dynamics on S is completely fixed and V/ G C°°(S) we have f = v i§^ ~ Qi§^- 
The energy (130) has the form 

E = \ (q\ + q 2 2 ) + ^i0i + ^202- (133) 

Instead of (131), (132) we can write 

f = {E,fl + (134) 

01 = 0, (135) 

02 = 0, (136) 

V/ G C°°(W). The reduction to S must be "canonical" in Hamiltonian formalism. The 
canonical transformation Pi = pi — q 2 /2, Q 1 = q 1 /2 + p 2 , P 2 = P2 — <7i/2 = 02, Q2 = 
g 2 /2 + Pl = 0i gives 5 = {P 2 = 0,g 2 = 0} and ft = P \ s = § (P? + Qf). The evolution 
equation (134) on S becomes / = {h, /}, / G C°°(S). 
Example 4- Now consider the Lagrangian (Ref . |]T6| ) 

21 



As it was in the previous example, W = TR 2 ($)T*R 2 , Ker uj = span {d/dvi,d/dv 2 }. The 
energy 

E = p 1 v 1 +p 2 V2-^e qi vl (137) 
produces (by (94)) the evolution vector field 

X E = + v 2 ^- + \e q2 vl-^- + + v 2 ^- (138) 

oqi oq 2 2 op 2 ov\ ov 2 

and the primary constrained form 

uj e = 4>\dvi + 4> 2 dv 2 = 0, (139) 

where 0i — P\ — e q2 Vi, 2 = p 2 . The condition (101) gives Vi = —V\V 2 as well as the 
new constraint <\$ — V\ — 0. Thus the primary constraint 0i is of the second class. On 
Mi = {x e W : 0i = 0} we have 

• df df 1 q2 2 d df df 

f = V!— + v 2 — + -e 92 Wi- v^ — + ^2^—, (140) 

oqi dq 2 2 ap 2 <7t>i dv 2 

02=^2 = 0, (141) 

0^ = «i = 0. (142) 

Since ~ the acceleration v 2 remains undetermined. The constraints 02 and 2 1 ' ) are 
of the first class. We have also 



E \ Ml = \p\e~ q2 + p 2 v 2 = h + v 2 (f) 2 , (143) 



and the system (140) -(142) in Hamiltonian form reads 



f = {hJ} + v 2 {<f )2 J} + v 2 ^ (144) 

02 = P2 = 0, (145) 
0^=^ = 0. (146) 



Once again, the acceleration v 2 remains arbitrary. 
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Example 5: Finally, consider the Lagrangian (Ref. Wh ) 

1 = 7^1 + q 2 Vi + (1 - a)qiv 2 + ^(qi - q 2 f 

The form u = dpi A dqi is degenerate on W = TR 2 ($T*R 2 , Ker u> = span {d/dvi,d/dv 2 }. 
The energy 

E = piVi + p 2 v 2 - I (147) 

gives the evolution vector field 

d d d d d d 

X E = V!— + v 2 — + (f 2 (l - a) +(3(q 1 - q 2 )) h («i - (3(qi - q 2 )) ^ Vv x - Vv 2 ^— 

oqi oq 2 dpi op 2 ovi ov 2 

(148) 

as well as the primary constraint form 

uj e = (ftidvx + (f) 2 dv 2} (149) 

where <p\ = p\ — q 2 — V\, <p 2 = p 2 — (1 — a)q\. For simplicity we consider here only some 
particular cases. 

Case A: a = 0,/3 = 0. 

The condition (101) gives v\ = 0, so (f)\ is the primary constraint of the second class. On 
M x = {x G W : 0i = 0} we get 

d d d . d . 

X El =vi- Vv 2 - Vv x - Vv 2 — . (150) 

c% dq 2 dp 2 dv 2 

The acceleration v 2 cannot be determined. The dynamics is controlled by 

/ = X El f, (151) 
02 = p 2 ~ qi = 0. (152) 



Since 



1 2 

E |mi= (vi + q 2 ) vi +p 2 v 2 - I = - (p! - g 2 ) + « 202 = + v 2 (p 2 , (153) 
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the Hamiltonian analog of (151) -(152) is 

f = {hJ} + v 2 {<j> 2 J} + v 2 ^-, (154) 
02 = P2 ~ 9i = 0. (155) 

Case B: a = 0,(3 ^ 0. 

In this case the compatibility condition allows to determine v\ = j3{q\ —q 2 ) and gives the 

secondary constraint 2 ^ = 92 — 9i- The primary constraint 0i is of the second class. On 

Ml = {x E W : 0i = 0} the evolution is given by 

9/ 9/ df df df 

f = Vitt- + v 2jr - + (ui - /3(9i - 9 2 )) 77— + P(qi - 92)77— + ^277—, (156) 
dqi dq 2 dp 2 dv 1 dv 2 

02 = P2 - 9i = 0, (157) 

0^ = 92 - 9i = 0. (158) 

The next iteration produces 2 — v 2 — V\ and the next one gives v 2 = (3{q\ — q 2 ). As a 
result, on S 3 C Mi, S 3 = {1 6 Mi : ^2 = 0,0^ = 0, 2 2 ) = 0}, we have / = and 
E \s 3 = \v\. To work out the Hamiltonian scheme, note that 



E |mx= ^ (Pi - 92) 2 - ^ (9l - 92) 2 + ^202 = /i + ^202- (159) 

The equation (145) then reads 

f = {hJ} + v 2 {<j> 2 J} + v 2 ^-. (160) 

Again, cf)^ = v 2 — v i = 0^ and 2 2 ^ = -u 2 = 0. Since t> 2 = — {/i, 2 L ' ) }, we can write (149) as 

/ = {h, /} - {/i, 0«} {0 2 , /} + {/i, 2 } {0«, /} (161) 

(we have added the last term "by hand" to antisymmetrize the bracket, since {h, 2 } ~ 

02^ = on the surface of the second class constraints. The canonical transformation of the 

form Pi = p 1 + p 2 - q 1 - q 2 , Q l = q u P 2 = p 2 - q l = 2 , Q 2 = q 2 - qi = 2 1) allows to 

(2) 

represent the dynamics on S3 = {P 2 = 0, Q 2 = 0, 2 = 0} in a simple way, / = {h, f}p 1) Q 1 , 
where h = Pf/2. 
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APPENDIX: THE GLOSSARY 

For convenience of the reader we give some notations and definitions of geometric objects 
used throughout the paper. More detailed information can be found in Ref. [11], [23| . 



1. M (or Q) denotes finite-dimensional real smooth manifold, (TM,tti, M) - its tangent 
bundle with local coordinates {q\v 1 }, (T*M, 7r 2 , M) - its cotangent bundle with local 
coordinates {q\Pi}, i = I, . . . N =dim M. 



2. X(M) denotes the Lie algebra of smooth vector fields on M. Locally X G X(M) is 



X = X*(q)£; } X*(q)eC°°(M). 



3. Correspondingly, Y G X(TM) in local coordinates reads Y = Y r (q, v)-^j + V' l (q, v)-^, 
where Y % and V % are smooth functions on TM. 

4. X V (TM)={X G X{TM) : vr,X = 0} is the algebra of vertical vector fields on TM. 
Any X G X V (TM) can be locally represented as X = X r (q, v)J~. 

5. Lift T X{M) -> X V (TM) maps X = X\q)^ G X(M) to X v = X\q)^ G 
X V (TM), where X v is the generator of the one-parameter group of the fiber dif- 
feomorphisms t — > (q l , v % + X l (q)i). 

6. The canonical (1,1) tensor r on TM is defined as the composition r = 7 o 7r*. In 



coordinates we have r = dq l * 



7. Cartan - Liouville one-form is defined for / G C°°(TM) by 9f — df o j. If — 9 
has maximal rank, cjj = d9f defines the symplectic structure on TM. 
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8. Liouville vector field A e X V (TM) is the generator of the one-parameter group of 
dilations (q,v) — > (q,e t v). Locally, A = f 1 ^-. 

9. Whitney sum E 1 E 2 of two bundles E 1 and E 2 with the same base B is the bundle 
with the base B and the fiber E lx (&E 2x Vrc e B. 
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